Studies of di-jet survival and surface emission bias in Au+Au collisions via angular 
correlations with respect to back-to-back leading hadrons 
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We report first results from an analysis based on a new multi-hadron correlation technique, ex- 
ploring jet-medium interactions and di-jet surface emission bias at RHIC. Pairs of back-to-back 
high transverse momentum hadrons are used for triggers to study associated hadron distributions. 
In contrast with two- and three-particle correlations with a single trigger with similar kinematic 
selections, the associated hadron distribution of both trigger sides reveals no modification in either 
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relative pseudo-rapidity, A77, or relative azimuthal angle, A<j>, from d+Au to central Au+Au col- 
lisions. We determine associated hadron yields and spectra as well as production rates for such 
correlated back-to-back triggers to gain additional insights on medium properties. 

PACS numbers: 25.75.-q 



Di-hadron correlation measurements in heavy ion colli- 
sions have advanced the studies of hot and strongly inter- 
acting matter at RHIC. The effects of jet-medium inter- 
actions are reflected in the broad away-side distributions 
of associated hadrons in di-hadron azimuthal correlations 
and their softened transverse momentum (pr) spectra 
4]. For associated hadrons with transverse momentum 
between 1.0 < px < 2.0 GeV/c in addition to the over- 
all broadening of the away-side {A<j> ~ tt) a dip at tt is 
reported in jj-ly]. Various scenarios and different mech- 
anisms of jet-medium interactions have been proposed 
to explain the experimental observations On the 

near-side (Acf> ~ 0) the modification of the correlation 
structure is also observed 0, [IlT - fl3l |. In [l2[ the same 
side peak is described by two components: the ridge, a 
long-range A77 plateau, and the small-A?y jet-like peak. 
It is argued that the ridge might arise from jet-medium 
interactions. However, recent three-particle analyses re- 
port possible decoupling of the jet and ridge features 14 1 . 
The existing three-particle correlation analyses do not 
specifically identify the away side jet hence the pseudo- 
rapidity correlation relative to the away-side jet cannot 
be studied. 

Di-hadron azimuthal correlations of two high-pr 
hadrons have been observed to exhibit jet-like peaks in 
both near- and away-sides [lj| [l6[ • Little modification 
is apparent in the shape of the peaks, but the ampli- 
tude of the away-side peak shows a strong suppression 
in central Au+Au collisions with respect to d+Au col- 
lisions. The data may be interpreted in two scenar- 
ios: in- medium energy loss followed by in- vacuum frag- 
mentation, and finite probability to escape the medium 
without interactions. As jet-medium interactions man- 
ifest themselves via shape changes of correlation struc- 
tures predominantly in the softer momentum range (be- 
low 2 GeV/c), we investigate the modification dynam- 
ics using a new 3-particle (2 + 1) correlation technique, 
which measures the degree of correlation of a softer par- 
ticle with a di-hadron trigger. The technique utilizes a 
second high-pr trigger on the away-side of the highest-pr 
particle in the event as an estimate of the away-side jet 
axis direction. 

The analysis presented in this letter used data col- 
lected by the STAR (Solenoidal Tracker at RHIC) ex- 
periment [T3| in the years 2003-2004. Data samples con- 
sist of 5.5M d+Au, 14M minimum bias Au+Au and 
18M central triggered Au+Au events from collisions with 
the center-of-mass per nucleon pair energy of 200 GeV. 
Based on the charged track multiplicity at mid-rapidity 
(\rj\ < 0.5) recorded by the STAR Time Projection 



Chamber (TPC) [18| the minimum bias Au+Au data 
are divided into 4 centrality bins. They correspond to 0- 
20%, 20-40%, 40-60%, and 60-80% of the total geometric 
cross-section. The central-trigger data represent approxi- 
mately 12% of the total cross-section. To express central- 
ity dependences we calculate via Monte-Carlo Glauber 
model the number of participating nucleons (^V par t) and 
the number of binary collisions (iVcou) for each central- 
ity bin used in the analysis. The tracks selection for the 
analysis is identical to [5j. 

We select the highest-pr particle in an event as the 
main (first) trigger and require Py Igl to be between 5 to 
10 GeV/c. The second trigger pr, p^ Ig2 , is selected to 
be above 4 GeV/c. We consider two high-px particles to 
be a di-jet-like pair if they conform to the requirement of 
relative azimuthal difference of |0 t rigi — 0trig2 — vr| < a. 
The a is chosen to be 0.2 corresponding approximately 
to lcr of the away-side peak width in relative azimuth 
for two high-pr hadrons [jjj . We construct the associ- 
ated hadron distributions in A(j> and A?/ coordinates with 
respect to each trigger direction. The associated soft par- 
ticle yield is measured relative to both triggers using the 
following expression 
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where iVt r jg is the number of trigger pairs, corrected 
for a dilution effect due to randomly associated pairs, 
and d 2 N ISiVJ / dArj dAtp is the associated hadron distribu- 
tion for all hadrons in triggered events. The associated 
hadrons p T (1.5 GeV/c < p^ ssoc < p^ rigl ) is selected to 
coincide with the range where the away-side dip in di- 
hadron correlations is reported. The "raw" distribution 
is corrected for single-track efficiency and acceptance ef- 
fects. The correction factor (e pa j r ) includes single track 
reconstruction efficiency, which depends on 77, <fi, pt and 
multiplicity. It also corrects for the pair acceptance ob- 
tained by the mixed-event technique. gAnSZwi re P re ~ 
sents the total background, consisting of multiple terms. 
The dominant background contribution is from combina- 
torics, and is accounted for by mixing trigger pairs and 
associated hadrons from different events, selected to have 
similar multiplicity and primary vertex position. The 
initially uniform distribution (apart from efficiency and 
acceptance effects) is then modulated by the flow term 
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where 
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II 



are sample-averaged values of 



the elliptical flow coefficient for primary and secondary 
triggers, and associated hadrons, respectively, to account 
for elliptical flow effects, which we consider independent 
of rj in the range studied 1|| 2^] ■ The multiplicity and px 



dependence of the elliptical flow coefficients (i^) for trig- 
gers and associated particles are obtained from averaging 
the results from the Event Plane and Four-particle Cu- 
mulant methods [l9j . The background due to randomly 
associated triggers in the initial selection of the trigger 
pairs is also considered @ . We use di-hadron correlations 
for the appropriate kinematic selections to determine the 
shape of the related background distribution, and the 
signal-to-noise ratio from the trigger-trigger correlation 
to estimate the per-trigger contribution for such pairs. 
The overall background level a zyam is estimated with the 



Zero- Yield at Minimum (ZYAM) method @, [H SI] . For 
the 2+1 correlation we choose the zero-yield region to be 
more than 3er (corresponding to f .3 radians) away from 
both jet-like peaks. We also use a double Gaussian plus a 
i>2-modulated background fit to estimate the systematic 
error due to the ZYAM method. To obtain the transverse 
momentum spectra for associated hadrons we follow the 
same procedure selecting the hadrons within 0.5 radians 
in relative azimuth and 0.5 in relative pseudorapidity of 
the respective trigger direction. 

Various sources of systematic uncertainties on mea- 
sured yields are evaluated. Single track reconstruction ef- 
ficiency contributes 5% overall normalization uncertainty 
for the associated hadron yields. This error does not 
affect same-side to away-side comparison for the corre- 
lated yields of the same data sample. The systematic 
uncertainty in the pair-acceptance correction due to fi- 
nite statistics of the mixed-events is estimated to be less 
than 5%. The uncertainty due to the elliptical flow cor- 
rection is estimated from the difference in the i>2 results 
from the Event-Plane and Cumulant methods. The un- 
certainty is found to be 5% in most central, less than 
1% in most peripheral and 9% in mid-central Au+Au 
events. This uncertainty is largely correlated between 
same- and away-side distributions, and is not applicable 
to d+Au events. Systematic uncertainty due to ZYAM 
normalization of the background level was estimated by 
varying the A<f> range for minimum yield region, and is 
found to be ~10% in the Au+Au data and less than 5% 
in d+Au events. The related uncertainty is correlated 
between same- and away-sides. Systematic uncertainty 
due to correlated background subtraction from the di-jet 
sample is determined to be less than 3% for both d+Au 
and Au+Au events. This error is correlated between 
the same- and away-sides, and was estimated by varying 
background normalization for the trigger-trigger correla- 
tion in a similar manner as for 2 + 1 correlation. Finally, 
the uncertainty in the di-hadron distributions used for 
estimating the correlated background term arises mainly 
from the uncertainty in and ZYAM normalization. We 



evaluate these in a similar manner and estimate the re- 
sulting effect to vary from less than 1% in d+Au events 
to about 5% in central Au+Au data. 

Figure 1 shows A</> and Ar\ projections (symmetrized 
about 0) of the correlation function defined above, mea- 
sured for central and mid-central Au+Au collisions and 
compared with the corresponding measurement from 
minimum bias d+Au data. The Ac/) projection (Fig. la) 
reveals jet-like peaks on both the same- and away-sides. 
The observed structures from Au+Au data of all central- 
ities are consistent within errors with the ones in d+Au 
collisions. This constitutes the first observation of the 
away-side correlation structure peaked in both A<^> and 
A?7 for the associated hadrons in this kinematic range. 
These results differ significantly from earlier di-hadron 
correlation observations, where substantial differences in 
both yields and shapes on both the same and away-sides 
were observed between Au+Au and d+Au data [5j. We 
observe no evidence of the Mach-cone effect in the three 
particle (2+1) correlations studied. We attempt to ex- 
plore the ridge properties by looking at A77 projections 
for di-jet triggered correlations for same-side as well as 
away-side triggers (Fig. 16). We find that the same-side 
integrated yields at 0.5 < Ar/ < 1.5 are consistent with 
zero; however, statistical limitations ban exclusion of the 
ridge magnitudes reported in [l2| . The 2 + 1 results pro- 
vide no evidence of a long-range "ridge" associated with 
the direction of either of the high-py triggers, consistent 
with lack of medium-induced effects on those di-jets se- 
lected by this analysis. The possible effects of radial flow 
on Mach-cone-like and rid ge s ignals in di-hadron correla- 
tions have been studied in [23J . We expect no such radial 
flow-related effects for tangentially emitted trigger pairs. 

We further explore the associated hadron spectra to 
complement the correlation shape analysis. Figure lc 
shows same- and away-side associated hadron yields plot- 
ted vs transverse momentum in d+Au and 12% central 
Au+Au collisions. Within the uncertainties no differ- 
ence is observed between associated hadron spectra on 
same-side or away-side, Au+Au or d+Au data samples. 
Earlier measurements from di-hadron correlations in a 
similar kinematic range have shown significant softening 
of the away-side associated hadron spectra |2J . 

To check for the possible away-side softening which 
is indicative of energy deposition in the medium, we 
have estimated the jet energy difference between same- 
and away- sides. We sum the transverse momentum 
of trigger particle and the charged hadrons used for 
the spectra as an estimate for jet energy in the fidu- 
cial range. We then calculate the difference between 
such energy estimates for the same-side and away-side 
for each of the data samples. The direct subtraction of 
the momentum- weighted same-side and away-side distri- 
butions used for this calculation results in a major im- 
provement of systematic uncertainties as no background 
re-normalization/subtraction is needed and uncertainties 
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FIG. 1: Projections of 2+1 correlation on A0 (o) and Arj (b) for 200 GeV top 12% central-triggered (squares) and mid-central 
(band) Au+Au and minimum bias d+Au (circles) data. Errors shown are statistical. The kinematic selection is as follows: 
5 < p^ lgl < 10 GeV/c, 4 GeV/c < p^ lg < p^ lgl , 1.5 < p^ soc < p^ lgl . c) Transverse momentum distributions per trigger pair 
for the same- and away-side hadrons associated with di-jet triggers (\Acj>\ < 0.5, |Atj| < 0.5). Errors are the quadrature sum of 
the statistical and systematic uncertainties. Inclusive charged hadron distribution 10% most central Au+Au data [24| is shown 
for comparison. 



on elliptical now cancel. We find that in the 12% central 
Au+Au data A(Au+Au) = S(p T ) samo - E(p T ) awa y = 
1.59 ± 0.19 GeV/c, similar to the minimum bias d+Au 
data value of A(d + Au) = 1.65 ± 0.39 GeV/c. The 
initial state kinematic effects are expected to cause simi- 
lar differences [3, [2(| , and calculated theoretically to be 
1.6 GeV/c. 

The absence of a jet quenching signal or a suppres- 
sion of associated hadrons suggests that the 2 + 1 parti- 
cle triplets considered in this analysis are biased towards 
surface jet emission, i.e. di-jet events where both jets are 
emitted nearly tangentially to the medium surface. In 
such a scenario, one would expect that di-jet production 
rates are determined by the surface area of the fireball. 
We investigate this possibility by studying the centrality 
dependence of jet and di-jet production rates in 200 GeV 
Au+Au data. 

Figure 2a shows the centrality dependence of the nu- 
clear modification factors u (ratio of binary-scaled 
per-event trigger counts in Au+Au and d+Au data) for 
the primary (single) triggers and di-jet triggers. The ob- 
served -Rj+Au U f° r single triggers is consistent with the 
suppression factors observed in inclusive charged hadron 
measurements |24j |. while correlated trigger pair produc- 
tion rates are suppressed even further. 

To examine if these results are consistent with the 
purely tangential jet production scenario, Monte Carlo 
Glauber calculations to find the spatial distribution of 
hard scattering were performed for each centrality bin. 
We assume a simplistic scenario where the interaction 
zone in heavy-ion collisions consists of a completely 
opaque core (full jet attenuation) surrounded by a per- 
meable corona (no jet-medium interactions). Energy loss 
fluctuations as discussed in I27| have not been taken into 
account. Participant eccentricity for the medium cores 



shape has been used. This has been found most suitable 
for elliptical flow results [28j |. 

The size of the core is estimated directly from the 
Raa suppression measurements 24, [29| and the calcu- 



lated core eccentricity. The total size of the interaction 
area is calculated in a similar way, requiring the inte- 
grated collision density inside this area to be 99% of the 
total N co i\. We further assume that all observed trig- 
ger hadrons come from the corona. Core/corona sizes 
that accommodate the inclusive hadron suppression re- 
sults are shown in Fig. 2b. The gray band shows propa- 
gation of uncertainties in the published Raa values into 
our calculations. Figure 2c displays the conditional di-jet 
survival rates. The expected rates are shown as a band, 
where the width reflects the uncertainty in the published 
Raa data. 

To estimate the conditional di-jet survival probabil- 
ity in the data, we use the double ratio I 



Au+Au 
d+Au 



Raa (trigger pairs) which rcflects changes in probabil- 

Haa (single triggers) ' J ° L 

ity to find an away-side trigger for each primary trigger 
in Au+Au data relative to d+Au. For 12% most central 
Au+Au data the likeliness of finding an away-side trigger 
for each primary trigger, which might be interpreted as 
conditional di-jet survival probability, is approximately 
20 % ±2(stat) ± 4.5(syst) compared to d+Au. This ra- 
tio is shown as symbols in Fig. 2c for comparison with 
the above model. However, we must point out that core 
emission where neither of the di-jets interacted with the 
medium cannot be ruled out by this analysis. 

To summarize: we attempt to investigate mechanisms 
of jet-medium interactions in ultra- relativistic heavy ion 
collisions using a novel technique involving three particle 
(2 + 1) correlations of hadrons associated with a corre- 
lated pair of back-to-back high px particles. Both same- 
side and away-side peaks of di-jet triggered correlations 
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FIG. 2: a) Relative production rates for jets and di-jets in Au+Au collisions of different centralities with respect to d+Au 
data, b) Calculations of the Glauber-based core/corona model that accommodates inclusive hadron suppression and single jet 
production rate results, c) Conditional di-jet survival probability in Au+Au data compared to d+Au reference. The band 
shows the expectations from Glauber-based core/corona model described in the text. Error bars in a) and c) are the quadrature 
sum of the statistical and systematic uncertainties. 



from central Au+Au data are found similar to the struc- 
tures observed in d+Au collisions at the same energy. We 
observe the same result in all centrality bins of minimum 
bias Au+Au data. We measure associated hadron spec- 
tra on each side of a di-jet primary trigger, and find that 
spectral shapes and integrated yields in central Au+Au 
data are similar on the same- and away -sides, and are 
consistent with those measured in d+Au data. In con- 
trast with earlier measurements with a single trigger, 
we observe no modification in correlations of associated 
hadrons in the momentum range of p^ ssoc > 1.5 GeV/c 
with the kinematic selection studied. Systematic assess- 
ment of di-jet production rates supports tangential emis- 
sion bias in a simplistic core/corona scenario, showing no 
evidence for di-jets interacting with the core. 
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